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Summary 


This report presents maps and tables for the global ocean tides, 69° N to 
68°S, derived from two years of Geosat altimetry. It supplements a paper 
by Cartwright and Ray which examines the energetics implied by the same 
data. Global maps of local and Greenwich admittance of the (altimetric) 
ocean tide, and maps of amplitude and Greenwich phase lag of the ocean tide 
are shown for M 2 , S 2 , N 2 , Oj, and Kj. Larger scale maps of amplitude and 
phase are also shown for regional areas of special interest. Spherical harmonic 
coefficients of the ocean tide through degree and order S are tabulated for 
the six major constituents. 
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1 Introduction 


The oceanic tide maps and spherical harmonic coefficients depicted or listed 
in this report were all derived from an improved (1990) analysis of two years 
of Geosat altimeter data, using techniques that are basically described in 
Cartwright and Ray (1990), a paper whose work was completed early 1989. 
The material herein forms a complement to a recent paper by the same 
authors ( Cartwright and Ray, 1991), which examines the energetics of the 
ocean tides implicit in the maps and compares dissipation rates with those 
derived from recent gravimetric analyses of satellite orbits, such as Marsh et 
al. (1990). The two papers by Cartwright and Ray will be referred to here 
as CR90 and CR91, respectively. 

CR90 was based on an analysis of the first year (1986-1987) of Geosat’s 
Exact Repeat Mission (ERM), and was restricted to latitudes between 60°N 
and 60°S. CR91 and the present compilation result from an improved pro- 
cessing of about twice the amount of ERM data and include extensions into 
the following sea areas which were not well covered, if at all, in CR90: 

(a) high latitude oceans between 68°N and S, including most parts of 
the Southern Ocean and of the Bering, Labrador, and Greenland- 
Norwegian Seas not permanently covered by ice; 

(b) marginal tidal seas omitted from or only partially represented in CR90, 
including the Mediterranean, Caribbean, Okhotsk, and Japan Seas, the 
South and East China Seas, and the broader part of the North Sea. 

2 Improved Data Analysis 

Our basic analysis of the altimetric sea surface with retention of basin-scale 
wavelengths is similar to that described in CR90. Here we give only a brief 
summary. We work exclusively with the variable 


a ap^) = ap„i)-a^t + r), o) 

where ( is the estimate at time t of the sea surface elevation at a dense 
network of earth positions P,, and r feJ 34.1 days is twice the orbital repeat 
period of Geosat. The elevations C are computed from the altimetry in 
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the usual manner, except that the computed “ocean tide” correction is not 
subtracted. The altimetric signal together with the “body tide” of the solid 
earth are subtracted from the orbital radius, computed to high precision, 
with corrections for ionospheric and tropospheric variability and for biases 
due to sea-state and atmospheric loading. The mission requirement that 
the Geosat groundtrack repeat to within ±1 km ( Bom et al., 1987) ensures 
elimination in (1) of all geographically fixed variables, including principally 
the geoid ( Brenner et al., 1990). The constancy of t to about 1 s imparts 
a simple constant transformation to every harmonic component of the tide, 
easily compensated for by imposing the same transformation on the tide- 
potential functions with which A( is correlated. 

Each estimate of £ is associated on the data tape supplied by NOAA with 
a nominal “ocean tide” Cmodei, computed from the Naval Surface Weapons 
Center (NSWC) model due to Schwiderski (1983). Denoting by Ar/ the 
sequence of height-differences similar to (1) but with ( replaced by y = 
C — Cmodei, we approximate the major part of the residual orbital radius error 
by least-squares adjustment of the function 

E(t) = A(t) cosujt + B(t) sinod + C(t) (2) 

to long time-sequences of &r) covering many orbital revolutions. In (2), u 
is the orbital frequency, about 89.91 rad. day -1 , / is referred to an ascending 
node, and A,B,C are slowly varying, modulatory functions, expressed as 
low-order polynomials in t. Finally, 

AC'(P„0 = AC(P„ t)-E(t) (3) 

is analyzed for its intrinsic tidal signal in time and surface position. This 
signal (apart from uncorrected errors and oceanographic noise) is what is 
usually called the “ocean tide,” slightly modified by elastic yielding of the 
sea bed, that is, by the “load tide” of the solid earth. 

Our procedure for tidal analysis is unorthodox but accurate in its speci- 
fication of the tidal potential and of a generally smooth admittance at every 
location. The interested reader is referred to CR90 for details. In brief, all 
values of A ('(P,,t) are collected in bins, 1° in latitude by 360/244 = 1.4754° 
in longitude, the latter being the spacing of adjacent ERM ground-tracks. 
Within each bin the data are convoluted with an array of 12 “orthotide” func- 
tions ( Groves and Reynolds, 1975), 6 representing the diurnal part of the tide 
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potential, 6 representing the semidiurnal part. The resulting 12 coefficients 
define smooth complex admittances in frequency space to the two species 
of potential. Maps of real and imaginary parts of the admittance at spe- 
cific frequencies, after some spatial smoothing between adjacent bins, define 
the spatial distribution of amplitude and phase of each Darwinian harmonic 
constitutent through use of the tables of Cartwright and Edden (1973). 

For the present study, the following improvements to the above scheme 
were carried out: 

1. The orbital radii were re-computed by Bruce Haines (University of 
Colorado) and Ron Williamson (STX) using the GEM-T2 gravity field, 
giving a global rms radial orbit error of about 19cm ( Haines ct at., 
1990), much lower than many previous computations for the Geosat 
orbit. The 6-day arcs used defined limits of the continuous sequences 
of A( (equation 1) which varied from 0.5 to 4.5 d without discontinuity. 
Over such time limits the modulatory functions /l, B , C of our once-per- 
revolution correction E (equation 2) could be adequately parameterized 
by linear forms Ao + A\t, etc. 

2. The total time span of the data series in t , starting 1986 November 8.0 
was increased from 34 Id to 6S5d plus the additional 34 d required to 
complete the definition of A<,\ This span includes a few short lacunae 
totalling 16 d but covers more than twice the alias-period (317 d) of 
the leading tide harmonic Mj at each grid position. The scatter of 
all tide parameters about smooth variations in longitude were thereby 
demonstrably reduced by a factor of about 1.5. 

Figures 1 and 2 illustrate this reduction in scatter for the particular 
case of the quadrature component of the N-j tide along a 1° zonal strip. 
Random scatter is most notably reduced in the quiet Pacific sector, 
170° 270°, but also in the noisy region of the Agulhas current system, 
15° 30°. An obvious outlier in the western South Atlantic at 317° was 
replaced later by interpolation between adjacent data. 

3. A number of import ant altimeter corrections were improved, as follows: 

(a) We used an improved empirical formula for sea-state bias, involv- 
ing both wave* height and wind speed (Ray and hoblinsky, 1991). 
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Figure 1: Raw tidal estimates of the quadrature component (//sin G) of the 
N 2 tide, along a 1° global strip centered at latitude 37.5°S, based on only 
the first year of Geosat ERM data. Each circle represents a separate tidal 
determination in a particular 1° x 1.475° bin. The solid line represents the 
Schwiderski (1983) model. The large scatter near longitudes 15°-30° is pre- 
sumably induced by the high mesoscale variability associated with the Agul- 
has Retroflection; similarly, scatter near longitude 310° reflects the presence 
of high variability at the confluence of the Brazil and Falkland Currents. 
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Figure 2: Similar to Figure 1, except based on two years of Geosat ERM 
data. The scatter is markedly reduced by the additional data, including 
the high scatter still somewhat present near the Agulhas and Patagonian 
high-mesoscale areas. These raw tidal estimates are spatially smoothed by 
methods discussed in CR90 before they are combined into grids to produce 
the maps shown below, e.g. Figure 5.2. A visual smoothing of the circled 
data shows that our N 2 component is several cm higher than Schwiderski’s 
(1983) estimates along the Pacific east of New Zealand and several cm below 
west of Chile. It is also 3-5 cm below Schwiderski across nearly the breadth 
of the South Atlantic at this latitude. 
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(b) We consistently applied the “inverse barometer” correction, which 
had been omitted from some areas of CR90. This correction was 
based on the Fleet Numerical Oceanography Center (FNOC) at- 
mospheric surface pressure data, which are provided on global 
grids every 12 hours. The corresponding Nyquist period of 1 day 
implies that atmospheric tides are not directly incorporated into 
the inverted barometer correction. The derived S 2 ocean tide 
showed marked improvement after this correction, but other con- 
stituents seemed little affected. 

(c) The use of water- vapor corrections from the FNOC model, instead 
of from the climatological model derived from Nimbus-7 data, gave 
us more comprehensive data access in marginal seas (see Introduc- 
tion, item b). 

After our data processing was completed, new sources for pressure and 
water vapor fields have become available, and some are thought prefer- 
able to the FNOC fields ( Ray et a/., 1991). But all atmospheric correc- 
tions are somewhat controversial; the robustness of tidal analysis relies 
as always on low correlation of residual errors with the tide potential. 

4. In order to reduce loss of data due to ice cover at latitudes above 60°, 
the collecting bins there were widened in longitude by a factor of 2 
or 3. During each satellite pass over a bin we characterized the several 
pulses by a median-filter instead of a straight mean value, to minimize 
the occasional effect of outliers. These procedures, together with the 
longer duration, greatly increased the number of usable passes for tidal 
analysis (above a practical minimum of 20 passes), resulting in fewer 
gaps and spatially smoother tide parameters. 
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3 Tidal Admittances 


As in CR90 and CR91, and earlier work, we have found it convenient to or- 
ganize the analysis and mapping of the tides in terms of complex admittance 
functions of frequency, defined as follows: Let a harmonic constituent of the 
tide potential of species m = 1 (daily) or 2 (half daily) at a surface point of 
latitude t p (N positive), east longitude A be expressed as 

V{<p, A, t) = gH 0 a™P™(p) cos(A <7 + at + mA), (4) 

where A a is the appropriate combination of classical mean (astronomical) 
longitudes at an arbitrary time origin, P™{<p) is the associated Legendre 
function, and H a is one of a set of known amplitudes in meters associated 
with a normalizing factor a™. In the normalization used for the tables of H# 
by Cartwright and Edden (1973), we have 

a- = —^(5/2470, ^(5/967r), 

P 2 m = 3 cos tp sin p>, 3 cos 2 <p, 

for m = 1,2 respectively. A tidal variation at P of the same frequency, 
expressed as 

h c {t) = H a (A'(<t) cos(A„ + at) - K(<7)sin(i4„ + <r0] ( 5 ) 

is said to have an admittance Z(a) = X(a) + iY(a) to the tide potential at 
the (Greenwich) origin of longitude. Its harmonic amplitude is 

H = \Z{a)\ H a (6) 

and its Greenwich phase lag is 

G = mw — arg (Z(<r)), (7) 

where the rmr arises from the negative sign of a 2 . 

The above reference to Greenwich transits is as arbitrary as the Greenwich 
meridian itself. In some applications, such as in CR91 where we were mostly 
concerned with calculations of the work done by gravity on the tides, it is 
more convenient to refer phases to local transit at the meridian of P, longitude 
A', say. In order to effect this we ‘rotate’ the admittance Z(a) to 

Z\a) — Z(a) exp(— im\'), (8) 
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so that 


X f (<7 ) = X cos m\ f + FsinmA' 
Y\o) = FcosmA' -XsinmAL 


Replacing Z by Z\ the amplitude relationship (6) still holds while the Green- 
wich phase lag G is replaced by the local phase lag /c, as originally defined by 
Kelvin and Darwin — see, for example, Schureman (1940, §§ 144-146). Thus, 

k = G + m\ f — mx — arg(Z'). (9) 

We shall call Z l (o) the local admittance in distinction to the Greenwich 
admittance Z(o). Its real part X* represents that part of the ocean tide at P 
which is in phase with the potential V at P; its imaginary part Y* represents 
the part in quadrature with V at P. 

Finally, in this survey of basic principles we must distinguish between 
the ocean surface tide as deduced directly from altimetry (minus the solid 
‘body tide’), which we shall call the ‘altimetric tide/ and the ‘ocean tide’ 
relative to the local sea bed as measured by in situ instruments. The alti- 
metric tide consists of the sum of the ocean tide and the ‘load tide’ of the 
solid earth induced by the ocean loading. In oceanic locations the altimetric 
tide is about 89-94 percent of the ocean tide, depending on the wavelength 
scale considered, but exact relationship requires an explicit Green’s function 
or high-degree spherical harmonic calculation. If the ocean tide is defined 
globally, the load tide and hence the altimetric tide may be computed for 
both ocean and land areas by multiplying the spherical harmonics of the 
ocean tide by a set of Love numbers ( e,g Ray and Sanchez , 1989). The 
reverse procedure is not so straightforward, because we do not know d priori 
the altimetric tide over land. An iterative procedure is then necessary to 
derive the ocean tide, as discussed in an appendix to CR91. (The iterative 
procedure converges so rapidly that, in fact, only one iteration is required.) 
We have applied this procedure to our altimetric tide solutions for all major 
harmonic constituents. The results, therefore, pertain correctly to the true 
ocean tide. 

True ocean tide maps are required for calculations of energetics and com- 
parison with hydrodynamic models. Altimetric tide maps are required for 
corrections to satellite altimetry prior to general oceanographic analysis. 
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Load tide maps are required for corrections to the positions of tracking sta- 
tions and similar precise geodetic applications. Examples of all three types 
of map will be presented in the following pages, but maps of the true ocean 
tides will be most comprehensive. Digital files containing arrays of all tidal 
types are available on request. 

An algorithm has been coded to compute the full diurnal plus semidiurnal 
altimetric tide (embodying all constituents through some 60 spectral lines) 
for any specified ocean position and time. It is also available on request. 
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4 Global Maps 

On the global scale, Figures 3. 1-3.4 show four pairs of maps for the M 2 
constituent, describing: 

3.1 contours of the local admittances X' and Y' of the altimetric tide; 

3.2 contours of the Greenwich admittances X and Y of the altimetric tide; 

3.3 contours of amplitude H and Greenwich phase lag G of the ocean tide; 

3.4 contours of amplitude Hi and Greenwich phase lag G\ of the load tide. 

Figures 4 . 1 - 4.4 show similar information for S 2 ; Figures 5 . 1 - 5. 4 show N 2 ; 
Figures 6. 1 - 6.4 show Oj; and Figures 7 . 1 - 7.4 show K,. In addition, the 
color Plate 1 depicts the differences in the in-phase HI = HcosG and 
quadrature H 2 = Hs'mG components of our latest M 2 map (Figure 3.3) 
and the corresponding values from Schwiderski’s numerical model. Color 
Plate 2 depicts similar information for Oj. 

The local admittances X Y' were also shown in CR91, and are repeated 
here for direct comparison with the Greenwich admittances X, Y, whose maps 
are similar to those of CR90 from the earlier analysis. Readers may find it 
interesting to compare the two types of representation, neither of which is 
common in the literature. In both types, nodes (amphidromes) occur at the 
intersection of zero contours of the real and imaginary parts, antinodes occur 
in areas of maximum modulus. The relation to (altimetric) amplitude and 
phase lag are 

H = H (X 2 + V' 2 ) 1 / 2 = H (A" 2 -(- V' ,2 ) 1/2 , 

G = mit — arg(.Y + iY), k = mw — arg(.Y ; + iY ' ), 

where II is the normalized potential amplitude from Cartwright and Edden 
(1973) (0.632m for M 2 , 0.294 for S 2 , 0.121 for N 2 , 0.262 for O u 0.369 for 
Kj), and m = 1 or 2 is the tidal species number. Note that the altimetric 
amplitudes are on average some 7 percent less than the true ocean tide for 
features of typical horizontal scales. As discussed in CR91, — Y'P™(<p) is 
directly proportional to the rates «_>f working of tidal forces. 

The maps of amplitude and phase (//, G) correspond to conventional 
ocean tide maps as produced by hydrodynamic models, such as those of 
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Schwiderski (1983) and Parke and Hendershott (1972). Differences from 
ground-truth at 80 selected stations, of order 2— 4 cm, are discussed in CR91. 

The load tide maps ( H t ,Gi ) may be compared with similar maps com- 
puted by Parke and Hendershott (1980) from their own model, and with 
load tides maps computed from Schwiderski’s ocean tide models by Ray and 
Sanchez (1989) and by Francis and Mazzega (1990). The maps shown here 
were computed directly from the altimetric tides as described in the appendix 
to CR91. The methodology exploits “fast” numerical transforms to decom- 
pose the tide into a high degree and order series of spherical harmonics (in 
this case, to degree and order 122 — see Section 6). The loading coefficients 
h' n required for these computations were taken from Farrell (1972); Ray and 
Sanchez (1989) note that the load tides are rather insensitive to the choice 
of h' n coefficients. The differences between the load tide maps given here and 
those in Parke and Hendershott (1980) and Ray and Sanchez (1989) stem 
primarily from differences in the underlying ocean tide models. 

The color prints in Plate 1 are presented for comparison with similar 
representations in CR90 from our earlier (less accurate) model. Note that 
the CR90 model was strictly altimetric, without correction for the load tide, 
whereas the present maps properly compare the ocean tide version from CR91 
with Schwiderski’s ocean tide. The discrepancies for Oi (Plate 2) have not 
been shown before in any form. Both M 2 and Oi from the CR91 model have 
lower rms residuals than Schwiderski’s models in comparison with the CR91 
ground-truth stations. S 2 and Kj, on the other hand, had slightly greater 
rms residuals than Schwiderski, so no color plates are shown for these solar 
constituents. 
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5 Regional Maps 

Figure pairs 8 to 14 show ocean tide maps in (H, G) notation for certain 
regional seas which are not clearly represented on the global scale, usually 
because of high amplitudes and steep gradients. Most show M 2 only, being 
the constituent of greatest interest, but the Ki tide is also shown for the 
Australasian Seas and the Sea of Okhotsk. 

The resonating Patagonian Shelf Sea (Figures 8a, b) has some extremely 
large amplitudes near the coast and two nodal points round which the tides 
rotate contra solem. (A third node off the coast of south Brazil, also included, 
is part of the main South Atlantic tide system.) Note that the amplitude 
contours are in steps of 10 cm from 0 to 60 cm, then in 20 cm steps upwards. 
Because of the large horizontal gradients and the proximity of land, including 
the Falkland Islands, some features may be poorly resolved or smoothed by 
the sampling procedure. Cartwright and Ray (1989) studied the energetics 
of this region from an early version of the CR90 analysis. 

The M 2 tide in the Mediterranean Sea (Figure 9) is well known to have 
mostly low amplitude; it is close to the threshold of resolution against noise 
from altimetric and orbit errors. Parts of the Adriatic Sea and the gulf be- 
tween 0°E and Gibraltar Strait are too narrow for adequate data capture 
and are not present in our data arrays. Therefore we do not show the nodal 
system in the northern Adriatic or the known rise of amplitude towards the 
Atlantic tide regime at Gibraltar. The figures show only the 10 cm ampli- 
tude and 120° and 180° phase lag, with indications of the intermediate values 
between isopleths. The rough division into 3 main basins (east, central, and 
west) with characteristic amplitudes and phases, together with resonant shal- 
low gulfs of higher amplitude in the Gulf of Gabes east of Tunisia and in the 
northwestern Adriatic, accord qualitatively with coastal data and the dy- 
namic model of Dressier (1980). Dressler’s amplitudes are < 10 cm between 
longitudes 0° and 10°E, whereas our data there are consistently in the range 
11-18 cm, probably too large in comparison to coastal data. 

The M 2 tide in the South and East China Seas (Figure 10) has a great 
deal of tight structure, especially in the northern part, known as the Yellow 
Sea. (For early but detailed maps based largely on coastal gauges, see Defant 
(1961, pp. 418-426); recent hydrodynamic models are described by Ye and 
Robinson (1983) and by Choi (1980).) There are four contra solem nodes in 
the Yellow Sea, but smoothing and the proximity of land tend to obscure 
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detail here. Nevertheless, the broad tidal features of the area appear to be 
correctly represented by the altimetry. 

Our M 2 maps for the Australasian seas (Figure 12) show the area of 
steep gradient of amplitude between western Australia and Java and the 
extent and limitations of altimetric data-retrieval in the island-studded seas 
of Indonesia. The large area north of Papua-New Guinea without phase 
contours, continued to the north in Figure 10, has phase lag everywhere 
between 270° and 300°, as documented elsewhere (e.g. Schwiderski , 1983). 
A nodal point is indicated inside the Gulf of Carpentaria, but our contouring 
algorithm has difficulty with such close-knit features in phase, emphasizing 
one advantage of (X,Y) plotting over ( H,G ) plotting. A node appears in 
the same Gulf for K 1? plotted in Figure 11. We have added a map for Kj to 
this region because of its well-known interplay of daily and half-daily tides 
between the Pacific and Indian Oceans. For general descriptions of tidal 
knowledge of Indonesian seas, see Defant (1961, pp. 417-421) and Bogdanov 
and Nefedyev (1961, 62). 

Turning finally to the northern part of the Asian Pacific border seas, Fig- 
ures 13a, b and 14a, b show respectively Kx and M 2 contours in the Okhotsk 
Sea area. The Okhotsk Sea is known for its strong diurnal response (Figure 
14.1); it also has large semidiurnal response and a number of semidiurnal 
nodes near the coast (Defant, 1961, pp. 430-433). The nodal points are not 
well resolved by the altimetry (Figure 14). The tides in the Japan Sea have 
low amplitude, confirmed by both Figures 14 and 10; few, if any, hydrody- 
namic models are available for the Japan Sea. Defant (1961, pp. 426-430) 
gives a general account of mid-century knowledge. 
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Figure 9: M 2 tide amplitude (top) and phase lag (bottom) in the Mediter- 
ranean Sea. Only 10 cm amplitude contours are shown, with + indicating 
> 10 cm and — indicating < 10 cm. Phase contours mark approximate po- 
sition of 120° and 180° phase lags, with figures indicating range of values 
between contours and coastline. 
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Figure 11: Ki tide in the Australasian seas. Amplitudes (top) in cm; Green 
wich phase lags (bottom) in degrees. 







Figure 12: M 2 tide in the Australasian seas. Amplitudes (top) in cm; Green- 
wich phase lags (bottom) in degrees. 
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Figure 13a: Kx tide in the Okhotsk Sea and adjacent seas. 
Amplitudes in cm. 
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Figure 14a: M 2 tide in the Okhotsk Sea and adjacent seas. 
Amplitudes in cm. 
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Figure 14b: M 2 tide in the Okhotsk Sea and adjacent seas. 
Greenwich phase lags in degrees. 
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6 Expansions in Spherical Harmonics 

For a number of applications it is desirable to have spherical harmonic ex- 
pansion coefficients for the major ocean tide constituents. Tables 2-7 below 
list the lowest degree and order spherical harmonic coefficients for the con- 
stituents Ox, Pi, Ki, N 2 , M 2 , and S 2 , computed from our global arrays of 
admittances, with adjustment for the load tide as explained in the Appendix 
to CR91 (i.e., the coefficients represent the pure ocean tide). These har- 
monics were computed through degree and order 122, the highest number 
appropriate for Geosat’s longitude spacing of 360°/244, but for space lim- 
itations, only those coefficients through degree and order 8 are presented 
here. 

As in Ray and Sanchez (1989), we have found it convenient to express 
the spherical harmonic series in the following convention: 

N N n 

\^2 a noPn( cos O) + cos m<f> - b nm smrn<j>)P™ (cos 0). (10) 

n=0 n=l m= 1 

The P™(fi) are the normalized associated Legendre functions, the first few 
of which are given in Table 1 . 

For any particular tidal constituent, the in-phase and quadrature compo- 
nents H cos G and H sin G are expanded as follows: 

N n , 

H cos G — (a nm cos m<f> — b nm sin m<f>)P™ (cos 0 ) 

n=0 m = 0 
N n f 

H sin G = ^ (c nm cos m<f> — d nm sin m<f>)P™(cos 0) 

n — 0 m=0 

where the prime on the second summation indicates that the first term is 
multiplied by y 2 , in agreement with (10). The coefficients a nm ,6 nm ,c„ m , and 
d nm are listed in the following tables in units of millimeters. 

Owing to conservation of mass, the terms a 0 o and Cq o should in theory be 
zero. All six listed constituents, however, are found to violate this condition, 
the most serious discrepancy being a qo for Kx. Nonconservation of mass 
in tide models is a common annoyance and it is here as likely induced by 
the neglect of seas in latitudes greater than 70° as by errors in the derived 
admittance maps. 
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Table 1: Normalized Associated Legendre Functions. 


^o°(z) = v/(l/2) 

/*(*) = v/(3/2) x 
F/(x) = v /(3/4)(l-x 2 ) 1 / 2 
j?(*) = v/( 5/2)(§x 2 -i) 

P 2 1 (x) = ^/l^SxCl-x 2 ) 1 / 2 

P 2 (x) = v/(5/48)3(l-x 2 ) 

^3°W = v/(7/2) (f* 3 -fx) 

Pi (*) = \/(7/24) (f x 2 - §)(1 - x 2 ) 1 / 2 
Pl( x ) = ^(7/240) 15x(l -x 2 ) 

P, 3 (x) = ^(7/1440) 15(1 - x 2 ) 3 / 2 


For comparisons to other studies, the following formulae convert our tab- 
ulated coefficients to the unnormalized coefficients often used in satellite 


geodesy: 


Dnm COS^ni m = \(“nm ± d nm )N™ 
Dnm sin V»nm = 

Dnm sin = — ^nm i d nm )N n 
Dnm COS V’nm = 2 ( C " m T 
where the normalization factor N™ is given by 

N m _ [ 2n + l(n-m)! ~ 
n 2 (n + m)! 


) 

I 


for species 2 


for species 1, 


The D± m and i/>± m are those used by Lambeck (1988), except that Lambeck 
adds 180° to ip for certain constituents like Oi and Pi whose amplitude in 
the potential expansion of Cartwright and Tayler (1971) is negative. Others 
have also used a phase angle £ — 90° — ip. Table 8 lists our coefficients in this 
D,ip convention for those terms that have recently been estimated through 
the precise tracking of artificial satellites ( c.g Marsh et al., 1990). 

In CR91 we discuss the harmonics of degree 2 in relation to their implied 
rates of tidal energy dissipation and in comparison with similar estimates 
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Table 2. Coefficients of Oi. 


Table 3. Coefficients of P x . 


n 

m 


^nm 

Cnm 

^nm 

0 

0 

6.423 


-3.685 


1 

0 

4.869 


-18.267 


1 

1 

1.621 

-7.434 

3.366 

-8.252 

2 

0 

8.691 


-12.353 


2 

1 

-15.747 

-41.610 

11.533 

-33.268 

2 

2 

16.599 

2.668 

-7.118 

17.479 

3 

0 

7.283 


15.837 


3 

1 

9.576 

16.476 

24.185 

31.418 

3 

2 

14.837 

30.911 

-11.574 

21.196 

3 

3 

-33.074 

-35.436 

2.886 

-8.137 

4 

0 

-3.856 


16.026 


4 

1 

-30.245 

9.224 

12.757 

-20.133 

4 

2 

6.781 

9.021 

-11.414 

8.549 

4 

3 

-6.373 

-7.109 

16.398 

-31.154 

4 

4 

7.852 

30.303 

1.954 

2.513 

5 

0 

-13.498 


5.191 


5 

1 

17.279 

-1.453 

-12.562 

12.811 

5 

2 

8.074 

2.626 

-6.827 

2.950 

5 

3 

5.629 

-7.539 

-0.215 

0.047 

5 

4 

5.679 

11.330 

-3.537 

4.161 

5 

5 

15.550 

5.789 

-0.855 

2.147 

6 

0 

-7.700 


4.749 


6 

1 

0.235 

-1.127 

-5.623 

0.233 

6 

2 

6.768 

-4.579 

0.574 

5.520 

6 

3 

-2.266 

-9.941 

1.193 

-1.077 

6 

4 

11.074 

7.313 

-5.935 

4.149 

6 

5 

0.960 

7.560 

2.361 

4.662 

6 

6 

-7.220 

-1.055 

-0.719 

-3.983 

7 

0 

-1.289 


-4.998 


7 

1 

0.758 

-0.679 

-10.129 

-4.609 

7 

2 

0.327 

-3.583 

5.312 

0.820 

7 

3 

-7.456 

2.589 

-3.849 

-5.772 

7 

4 

6.068 

7.049 

0.296 

7.041 

7 

5 

-5.400 

3.281 

-3.588 

-0.854 

7 

6 

-5.980 

-6.772 

1.614 

-5.693 

7 

7 

7.311 

-1.009 

-1.612 

0.123 

8 

0 

-1.499 


-8.442 


8 

1 

8.357 

4.993 

-2.736 

10.055 

8 

2 

-1.748 

-0.171 

3.987 

1.680 

8 

3 

-4.318 

-2.556 

8.497 

-2.909 

8 

4 

5.305 

8.476 

-1.392 

4.959 

8 

5 

-8.806 

-6.106 

0.597 

-4.211 

8 

6 

-8.487 

-1.447 

0.606 

-6.996 

8 

7 

0.432 

-1.851 

-0.913 

2.861 

8 

8 

3.191 

-8.136 

0.066 

3.617 


n 

m 

®nm 

^nm 

Cnm 

dnm 

0 

0 

4.339 


0.341 


1 

0 

2.813 


-9.486 


1 

1 

1.353 

-2.584 

2.962 

-4.268 

2 

0 

-2.729 


-4.435 


2 

1 

-4.522 

-14.256 

3.927 

-15.024 

2 

2 

10.181 

0.327 

-0.438 

12.517 

3 

0 

4.436 


10.126 


3 

1 

-6.182 

1.132 

13.637 

9.998 

3 

2 

2.725 

10.220 

-0.232 

17.350 

3 

3 

-13.200 

-14.972 

-2.517 

-11.932 

4 

0 

0.413 


3.876 


4 

1 

-13.111 

8.873 

0.281 

-7.227 

4 

2 

4.256 

1.322 

-5.597 

2.801 

4 

3 

-6.080 

-0.610 

5.812 

-11.886 

4 

4 

4.103 

17.736 

5.547 

5.212 

5 

0 

-8.082 


-0.946 


5 

1 

7.135 

-2.615 

-5.373 

5.182 

5 

2 

4.452 

-0.169 

-2.992 

2.334 

5 

3 

1.061 

-2.456 

0.201 

-1.580 

5 

4 

1.191 

2.701 

-3.836 

4.664 

5 

5 

7.059 

4.695 

1.152 

0.560 

6 

0 

-2.630 


3.039 


6 

1 

1.147 

-1.779 

-3.200 

1.223 

6 

2 

1.884 

-4.648 

0.577 

1.942 

6 

3 

-2.375 

-5.096 

-0.734 

-1.103 

6 

4 

6.361 

0.344 

0.519 

3.453 

6 

5 

-1.389 

3.758 

-1.366 

2.305 

6 

6 

-1.819 

0.579 

-4.790 

-2.057 

7 

0 

-0.035 


-0.822 


7 

1 

1.554 

0.871 

-4.209 

-1.893 

7 

2 

-1.395 

-1.943 

3.597 

-0.916 

7 

3 

-2.585 

0.985 

-1.583 

-4.014 

7 

4 

4.098 

0.995 

1.149 

5.058 

7 

5 

-3.028 

-0.578 

-1.894 

1.401 

7 

6 

-2.036 

-2.931 

-0.611 

-4.382 

7 

7 

4.359 

-2.228 

0.533 

0.840 

8 

0 

2.549 


-3.366 


8 

1 

4.339 

1.178 

-0.259 

3.721 

8 

2 

-1.485 

-0.402 

0.919 

0.595 

8 

3 

-3.329 

-0.271 

2.894 

-2.763 

8 

4 

2.646 

1.588 

0.952 

4.352 

8 

5 

-3.558 

-2.587 

-1.234 

-2.697 

8 

6 

-2.669 

-0.909 

0.239 

-1.161 

8 

7 

-0.019 

-1.919 

-0.805 

3.426 

8 

8 

1.011 

-4.174 

1.487 

0.653 


Units are millimeters. 
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Table 4. Coefficients of Kj . 


Table 5. Coefficients of N 2 . 


n 

m 

^nro 

bum 

Cfim 

dnm 

n 

771 

®nm 

^nm 

Cfim 

d-nm 

0 

0 

13.791 


0.353 


0 

0 

0.124 


-0.613 


1 

0 

8.866 


-28.802 


1 

0 

2.347 


-0.356 


1 

1 

2.508 

-7.660 

9.532 

-13.986 

1 

1 

-2.684 

0.900 

0.889 

-4.602 

2 

0 

-9.604 


-14.703 


2 

0 

-0.330 


-2.724 


2 

1 

-12.854 

-42.676 

11.208 

-45.238 

2 

1 

-5.427 

-0.315 

1.936 

0.744 

2 

2 

31.138 

-0.186 

-1.032 

38.567 

2 

2 

-9.838 

-17.890 

16.799 

-16.710 

3 

0 

13.331 


31.474 


3 

0 

1.807 


-11.885 


3 

1 

-21.039 

3.302 

39.674 

28.776 

3 

1 

7.518 

4.513 

11.610 

-7.367 

3 

2 

6.491 

28.839 

0.196 

53.202 

3 

2 

1.616 

6.610 

-4.544 

2.724 

3 

3 

-39.873 

-44.973 

-8.232 

-37.442 

3 

3 

25.920 

-3.790 

1.572 

1.172 

4 

0 

1.880 


12.224 


4 

0 

-12.528 


5.499 


4 

1 

-38.614 

27.833 

0.126 

-20.796 

4 

1 

4.211 

0.708 

10.177 

2.911 

4 

2 

13.531 

4.182 

-17.338 

7.734 

4 

2 

19.048 

25.642 

-6.994 

6.637 

4 

3 

-18.582 

-1.790 

16.866 

-35.992 

4 

3 

-16.758 

17.226 

-2.399 

10.398 

4 

4 

11.582 

54.126 

17.840 

17.097 

4 

4 

-21.380 

-12.878 

10.398 

-4.917 

5 

0 

-24.332 


-4.160 


5 

0 

-6.670 


10.245 


5 

1 

21.768 

-8.420 

-16.458 

15.317 

5 

1 

-18.060 

-12.981 

-13.043 

-5.285 

5 

2 

13.758 

-1.053 

-9.190 

6.908 

5 

2 

-18.941 

-7.079 

12.244 

20.830 

5 

3 

2.588 

-6.615 

0.373 

-5.212 

5 

3 

-11.707 

-9.867 

3.576 

-12.201 

5 

4 

3.881 

7.064 

-12.324 

14.223 

5 

4 

8.180 

-4.602 

8.199 

23.779 

5 

5 

21.345 

14.730 

4.217 

1.840 

5 

5 

12.910 

-13.687 

-6.123 

-2.481 

6 

0 

-8.082 


9.589 


6 

0 

14.584 


-6.514 


6 

1 

3.962 

-5.982 

-9.715 

3.940 

6 

1 

-2.806 

0.950 

-9.253 

-1.520 

6 

2 

5.456 

-14.469 

1.448 

5.199 

6 

2 

1.128 

-14.889 

5.353 

-2.409 

6 

3 

-7.520 

-15.453 

-2.455 

-3.920 

6 

3 

-1.457 

-0.653 

8.699 

0.641 

6 

4 

19.138 

0.145 

2.567 

10.136 

6 

4 

12.146 

2.412 

6.385 

1.645 

6 

5 

-3.804 

11.730 

-4.686 

6.979 

6 

5 

-3.179 

8.960 

-1.089 

-14.114 

6 

6 

-4.246 

1.946 

-15.247 

-6.344 

6 

6 

8.130 

14.598 

21.618 

8.964 

7 

0 

-0.415 


-1.870 


7 

0 

10.444 


1.792 


7 

1 

4.944 

2.691 

-12.436 

- 5.782 

7 

1 

3.632 

5.663 

-0.594 

4.968 

7 

2 

-4.688 

-5.724 

10.842 

-3.313 

7 

2 

3.422 

-5.350 

-8.075 

-13.139 

7 

3 

-7.873 

3.337 

-5.025 

-12.540 

7 

3 

-8.559 

-2.421 

-0.531 

-5.251 

7 

4 

12.493 

2.297 

3.860 

15.291 

7 

4 

-1.021 

-1.091 

-3.435 

8.477 

7 

5 

-9.324 

-2.758 

-5.880 

4.204 

7 

5 

-4.462 

5.581 

-5.225 

13.202 

7 

6 

-5.335 

-8.732 

-2.147 

-13.607 

7 

6 

-1.164 

-6.762 

-3.062 

6.700 

7 

7 

13.412 

-7.523 

2.131 

2.680 

7 

7 

9.327 

10.808 

2.891 

3.865 

8 

0 

8.315 


-9.542 


8 

0 

-1.330 


3.035 


8 

1 

13.244 

3.760 

-0.498 

11.048 

8 

1 

3.870 

1.355 

-4.851 

1.686 

8 

2 

-4.307 

-1.287 

2.521 

1.927 

8 

2 

-4.560 

2.444 

-4.959 

-1.868 

8 

3 

-9.822 

-0.300 

8.098 

-8.486 

8 

3 

6.601 

-4.516 

-3.579 

1.730 

8 

4 

7.669 

3.921 

3.186 

13.352 

8 

4 

-2.387 

-2.322 

11.113 

6.069 

8 

5 

-10.402 

-7.759 

-3.999 

-8.376 

8 

5 

-6.043 

6.228 

-3.443 

3.352 

8 

6 

-7.957 

-3.043 

0.710 

-3.175 

8 

6 

-4.416 

2.961 

-3.046 

2.010 

8 

7 

-0.006 

-6.373 

-2.448 

10.674 

8 

7 

-12.205 

5.693 

0.784 

-2.092 

8 

8 

2.567 

-12.637 

4.966 

1.528 

8 

8 

0.554 

-1.260 

-1.873 

-0.146 
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Table 6. Coefficients of M 2 . 


Table 7. Coefficients of S 2 . 


n 

m 

^nm 

bnm 

c nm 

^nm 

n 

m 


^nm 

c nm 

dnm 

0 

0 

6.520 


-2.747 


0 

0 

3.374 


3.269 


1 

0 

7.320 


5.956 


1 

0 

4.664 


2.571 


1 

1 

— 12.735 

15.518 

-9.016 

-12.699 

1 

1 

-2.498 

7.825 

-0.167 

0.292 

2 

0 

-7.225 


-30.398 


2 

0 

11.213 


-7.474 


2 

1 

-24.297 

6.794 

-6.219 

-11.706 

2 

1 

-5.644 

9.722 

-4.194 

1.099 

2 

2 

-56.604 

-73.104 

86.870 

-90.602 

2 

2 

-23.047 

-19.205 

38.661 

-21.652 

3 

0 

9.639 


-46.005 


3 

0 

6.892 


-26.029 


3 

1 

27.096 

23.199 

74.606 

-32.045 

3 

1 

-9.660 

9.082 

27.785 

-0.534 

3 

2 

-1.523 

23.062 

— 10.671 

9.774 

3 

2 

-5.958 

13.594 

-22.113 

-4.515 

3 

3 

108.486 

-22.920 

40.788 

-3.495 

3 

3 

21.705 

0.721 

26.562 

-14.133 

4 

0 

-40.916 


23.428 


4 

0 

-28.155 


2.737 


4 

1 

1.088 

-8.575 

63.547 

27.841 

4 

1 

-8.734 

-12.709 

26.185 

-5.349 

4 

2 

74.350 

94.796 

-18.889 

75.479 

4 

2 

37.375 

9.557 

-4.417 

31.725 

4 

3 

-70.380 

77.704 

-52.820 

71.706 

4 

3 

3.252 

7.504 

-20.393 

44.780 

4 

4 

-96.299 

-73.162 

11.401 

-37.282 

4 

4 

-20.990 

-29.871 

-7.898 

-30.279 

5 

0 

-19.777 


22.186 


5 

0 

-2.756 


21.254 


5 

1 

-69.804 

-37.849 

-100.333 

-45.503 

5 

1 

-1.100 

-0.981 

-47.657 

-16.770 

5 

2 

-81.484 

-56.592 

14.114 

86.186 

5 

2 

-17.797 

-25.752 

5.908 

30.786 

5 

3 

-40.991 

-25.305 

9.861 

-75.385 

5 

3 

-18.360 

14.777 

11.005 

-29.203 

5 

4 

27.995 

-35.100 

74.224 

97.607 

5 

4 

-17.051 

-12.240 

23.998 

50.654 

5 

5 

57.708 

-40.326 

-23.676 

-38.304 

5 

5 

28.951 

3.534 

-22.061 

-10.434 

6 

0 

56.985 


-9.897 


6 

0 

19.850 


-2.247 


6 

1 

-5.093 

-1.029 

-42.678 

-6.807 

6 

1 

5.445 

6.059 

-18.304 

2.912 

6 

2 

-15.925 

-65.591 

44.290 

-39.668 

6 

2 

-23.504 

-13.049 

5.107 

-27.751 

6 

3 

-20.135 

0.184 

46.227 

-3.366 

6 

3 

-19.845 

-0.550 

3.326 

11.428 

6 

4 

30.601 

-0.888 

36.936 

12.123 

6 

4 

8.190 

-4.832 

4.634 

-4.820 

6 

5 

-11.485 

54.338 

-6.474 

-55.435 

6 

5 

-9.713 

27.483 

8.027 

-20.275 

6 

6 

9.580 

56.721 

111.845 

78.788 

6 

6 

-19.454 

-4.257 

37.629 

44.844 

7 

0 

33.830 


33.420 


7 

0 

-2.808 


8.776 


7 

1 

20.385 

17.801 

10.123 

35.489 

7 

1 

2.416 

-2.663 

6.438 

10.629 

7 

2 

28.298 

-0.382 

-36.705 

-59.548 

7 

2 

13.437 

11.666 

-9.943 

-16.182 

7 

3 

-32.679 

-14.754 

-19.272 

-26.310 

7 

3 

-10.349 

-3.113 

-12.029 

-25.536 

7 

4 

-3.001 

-27.161 

-9.241 

36.102 

7 

4 

2.044 

-10.613 

-4.555 

2.185 

7 

5 

-12.942 

6.652 

-32.701 

74.351 

7 

5 

-0.261 

-8.221 

— 14.512 

24.055 

7 

6 

-7.486 

-42.845 

-8.875 

19.921 

7 

6 

-5.307 

-12.561 

-7.946 

-10.137 

7 

7 

48.794 

43.414 

31.749 

28.019 

7 

7 

8.356 

12.098 

22.316 

15.170 

8 

0 

-7.415 


12.428 


8 

0 

-0.156 


8.046 


8 

1 

31.814 

7.927 

-23.788 

5.829 

8 

1 

16.608 

-0.693 

0.276 

1.876 

8 

2 

-14.743 

15.690 

—36.526 

-14.510 

8 

2 

3.131 

6.705 

-13.693 

-6.271 

8 

3 

34.484 

-23.014 

— 14.665 

7.803 

8 

3 

19.426 

-10.977 

0.925 

-0.333 

8 

4 

-28.656 

-19.711 

43.145 

24.483 

8 

4 

-12.510 

-12.454 

7.358 

5.245 

8 

5 

-19.204 

22.334 

-28.312 

31.960 

8 

5 

0.515 

-1.337 

-14.901 

14.836 

8 

6 

-12.206 

5.368 

—23.358 

2.721 

8 

6 

1.686 

-0.975 

-8.086 

0.032 

8 

7 

-58.555 

21.136 

-17.010 

1.044 

8 

7 

-16.761 

0.108 

-14.250 

1.701 

8 

8 

19.809 

-5.813 

-10.962 

-7.817 

8 

8 

16.856 

1.118 

4.850 

-2.800 
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Table 8: Some selected harmonic coefficients. 


71 

Diurnals 
m D+ m (cm) 

tl& 

Tide 

n 

Semi diurnals 
™ Dtm (cm) 

< nm 

2 

1 

2.33 

42.7 

n 2 

2 

2 

0.70 

127.4 

3 

1 

1.13 

280.6 


3 

2 

0.10 

291.3 

4 

1 

1.20 

86.0 


4 

2 

0.23 

308.2 

5 

1 

0.69 

249.7 


5 

2 

0.08 

84.4 

6 

1 

0.09 

185.9 


6 

2 

0.06 

93.6 

7 

1 

0.19 

157.8 


7 

2 

0.03 

195.7 

8 

1 

0.34 

247.2 


8 

2 

0.02 

229.0 

2 

1 

0.86 

47.1 

m 2 

2 

2 

3.51 

132.6 

3 

1 

0.35 

343.0 


3 

2 

0.30 

283.7 

4 

1 

0.52 

112.9 


4 

2 

1.05 

322.8 

5 

1 

0.27 

257.4 


5 

2 

0.29 

86.2 

6 

1 

0.05 

239.0 


6 

2 

0.38 

116.8 

7 

1 

0.09 

176.2 


7 

2 

0.12 

229.3 

8 

1 

0.14 

259.9 


8 

2 

0.12 

240.7 

2 

1 

2.56 

47.2 

s 2 

2 

2 

1.18 

127.7 

3 

1 

1.00 

348.0 


3 

2 

0.32 

253.6 

4 

1 

1.55 

115.0 


4 

2 

0.39 

348.6 

5 

1 

0.81 

257.8 


5 

2 

0.14 

67.7 

6 

1 

0.17 

244.7 


6 

2 

0.17 

160.5 

7 

1 

0.28 

176.8 


7 

2 

0.05 

262.8 

8 

1 

0.42 

260.1 


8 

2 

0.04 

261.2 


from other models. As a minor detail, we point out in particular that the in- 
phase component of the (2,2) harmonic, namely .D^cos^k, ma y ' n error 
by 0.01 x the potential amplitude divided by g, on account of the subtraction 
of the body tide with an older version of the Love number h 2 (0.619 instead of 
the modern 0.609). The difference in M 2 caused by this correction is less than 
1 mm, significant in comparison with D 22 itself, but negligible in a general 
comparison of tidal amplitudes. We have not applied this correction to the 
values listed in Tables 2-8; however, the correction was applied to M 2 and 
S 2 in Table 5 of CR91. 
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Plate Captions 


Plate 1: Differences (cm) between the present Geosat-derived ocean tide 
and the Schwiderski (1983) model, for the M 2 constituent. (Top) In-phase 
component H cos G; (bottom) quadrature component HsinG. Each color 
spans an interval of 5 cm. Differences on the Patagonian Shelf exceed 50 cm. 


Plate 2: As for Plate 1, but for the Ox constituent. Each color spans an 
interval of 2 cm. 
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